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We study black hole singularities in the AdS/CFT correspondence. These singulari- 
ties show up in CFT in the behavior of finite-temperature correlation functions. We first 
establish a direct relation between space-like geodesies in the bulk and momentum space 
Wightman functions of CFT operators of large dimensions. This allows us to probe the 
regions inside the horizon and near the singularity using the CFT. Information about the 
black hole singularity is encoded in the exponential falloff of finite-temperature correla- 
tors at large imaginary frequency. We construct new gauge invariant observables whose 
divergences reflect the presence of the singularity. We also find a UV/UV connection that 
governs physics inside the horizon. Additionally, we comment on the possible resolution 
of the singularity. 
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1. Introduction 

The AdS/CFT correspondence [[j],|],|3[ provides exciting avenues for exploring various 
issues in quantum gravity. An important question that the AdS / CFT correspondence may 
shed light on is that of the nature of spacelike singularities, like the Big Bang or the singu- 
larity of a Schwarzschild black hole. A good laboratory for studying spacelike singularities 
is an eternal black hole in anti-de Sitter (AdS) spacetime. It has been conjectured that 
quantum gravity in an AdS^+i black hole background is described by a boundary confor- 
mal field theory on S d ~ 1 x IR at a temperature given by the Hawking temperature of the 
black hole @,|J|]0. 

The conventional wisdom regarding singularities is that they signal the breakdown 
of classical gravity and should go away when stringy or quantum gravitational effects are 
taken into account. Since in AdS/CFT, classical gravity corresponds to the large N and 
large 't Hooft coupling limit of the boundary theoryil, one expects that finite N or finite 
't Hooft coupling effects may resolve these singularities 0||. Such considerations suggest 
the following strategy: 

1. Identify manifestations of the black hole singularity in the large N and large 't Hooft 
coupling limit of the finite temperature boundary theory; 

2. From these manifestations, understand the precise physical mechanism through which 
the finite N or finite 't Hooft coupling effects may resolve the singularity. 

In the boundary theory, the physical observables are correlation functions of gauge in- 
variant operators. This means that the physics of singularities should be encoded in the 
behavior of boundary correlation functions in appropriate limits. 

One of the obstacles! in understanding black hole singularities from finite temperature 
boundary theory is that the singularities are hidden behind event horizons. The boundary 
conformal field theory evolves through the bulk Schwarzschild time, i.e. from the point 
of view of an external observer, and does not appear to directly describe the physics 
beyond the horizon. In other words, if time evolution inside the horizon of a black hole 
is to be described by the boundary theory, time has to be holographically generated. This 

1 More precisely, an AdS black hole appears as a saddle point in the path integral of boundary 
theories in the large N limit. 

2 While we are interested in black holes in all dimensions, when talking about boundary the- 
ories, we will more specifically have in mind d = 4, which are SU{N) Yang-Mills theories. 

3 in addition to the standard difficulties to decode the local bulk physics from boundary theories 
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makes the problem particularly challenging, while at the same time exciting. In particular, 
understanding physics beyond the horizon should shed light on how to holographically 
describe a Big Bang cosmological spacetimeS. 

A number of authors J7[|l7] have explored how to extract physics beyond the horizon 
from the boundary theory correlation functional. In particular, Fidkowski et al [12 1 found 
an interesting but subtle signal of the singularity in the boundary correlators. They found 
that AdSd+i black holes with dimension d > 3 contain spacelike geodesies, connecting two 
asymptotic boundaries, which could get arbitrarily close to the singularity. The authors 
further argued that such geodesies imply the presence of poles on secondary sheets of 
the analytically continued coordinate space correlation functions in the large operator 
dimension limit. 

Here we further explore the manifestations of the black hole singularity in the bound- 
ary theory and discuss their implications for the resolution of the singularity. 

We will establish a direct relation between space-like geodesies in the bulk spacetime 
and the large operator dimension limit of momentum space Wightman functions in the 
boundary theory. We show that physics in the region beyond the horizon is encoded in 
the behaviors of boundary correlation functions along the imaginary frequency axis. In 
particular, this gives a clear indication that the "time" inside the horizon is holographically 
generated from the boundary theory. The presence of the curvature singularity leads to 
certain exponential falloff of the correlation functions near the imaginary infinity. We also 
construct new gauge invariant observables which have singularities precisely reflecting the 
curvature singularity of the black hole. 

In this paper we will present the main idea using the example of an AdSs Schwarzschild 
black hole, leaving technical details and more extensive discussions of other examples to a 
longer companion paper ||23|| . In []23] we develop a new method for computing the large mass 
quasi-normal frequencies of the black hole, whose knowledge is crucial for the discussion 
of this paper and ]23| . 



The plan of the paper is as follows. In section 2, we review the relevant black hole 
geometry and the computation of boundary Wightman function in AdS/CFT. In section 



4 The regions around the past and future singularities of a black hole can also be viewed as 
Big Bang or Big Crunch cosmologies. 



See also [18,19|- For other recent discussion of spacelike singularities in AdS/CFT, see e.g. 

I • 
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3, we establish a connection between bulk spacelike geodesies and boundary Wightman 
functions in the large operator dimension limit. In section 4, we study the manifestation 
of the singularities in Yang-Mills theory. We conclude in section 5 with a discussion of the 
possible resolution of singularities at finite N. 



2. Boundary correlation functions from AdS/CFT 

2.1. Black hole geometry 

We will consider big black holes in AdSs, which have a positive specific heat and are 
the dominant contribution to the thermal canonical ensemble of the boundary Yang-Mills 



theory when the temperature is sufficiently high |25[ . 

The metric for a Schwarzschild black hole in an AdSs spacetime is given by 



ds 2 = -f(r)dt 2 + /(r)" 1 ^ 2 + r 2 dfl 2 3 



with 



f(r) 



r 2 + 1 



V- 1 / 2 

1 + 



r 2 )(r 2 + r 2 ), 



(2.1) 



(2.2) 



2 2 

r Q r x 



where \i is proportional to the mass of the black hole and the event horizon is at r = tq. 
ro, T\ can be solved in terms of fi. We have set the curvature radius of AdS to be unity, as 
we will do throughout the paper. As r — » oo, the metric goes over to that of global AdS 
with t identified as the boundary time. The fully extended black hole spacetime has two 
disconnected time-like boundaries, each of topology S 3 x M. 




Fig. 1: Penrose diagram for the AdS black hole. There are two asymptotic AdS 
regions, which are space-separated from each other. A null geodesic going from the 
boundary to the singularity is indicated in the figure. 
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It is often convenient to use the tortoise coordinate 



DC 



z= I #v = -#log(^)+A tan -i!l ( 2 .3) 



f(r) 4n \r + r J 2n r 
with 

jj = -£V (2.4) 

The region outside the horizon corresponds to z G (0, +oo) with z = at the boundary 
and z — ► +oo at the horizon. /3 is the inverse Hawking temperature. To have a feeling of 
the physical meaning of f3 we note that the complex Schwarzschild time that it takes for a 
radial null geodesic to go from the boundary to the singularity is given by 

± I~W) =± \^ ) (2 ' 5) 

where the imaginary part of the integral arises by going around the pole at r = r$ in the 
complex r-plane. A nonzero (3 implies that the Penrose diagram for the black hole is not 
a square, as was first pointed out in It is convenient to introduce a complex quantity 



B = P + l [3= 27r{r 2 1 + l 2 ro) =-gL. ( 2 . 6 ) 
Tq + r i n — itq 

which will be important in our discussion. One can invert ( f2.3| ) to find r(z). In particular, 
for Rez > 4, r is a one-to-one periodic function of z with period i^. 

We will consider the complexified Kruskal spacetime in which points related by 

■ m + n n .m — ti _ r , , n „ s 

t->t + i 0, z^z + i 0, m,neZ (2.7) 

are identified. The Lorentzian section (fig. |1|) of the complexified spacetime can be conve- 
niently described using (i, z) with constant imaginary parts. For example, up to identifi- 
cations (|2.7p , region III can be specified by 

Imi = -— , lmz = 0. (2.8) 

2.2. Boundary Wightman functions 

Now consider an operator O in the boundary theory corresponding to a bulk field <p 
of mass m. In the supergravity limit, the conformal dimension of O is given by [|3]|| 



±+v, u=^+m*. (2.9) 



Thermal boundary two-point functions of O can be obtained from free bulk Green functions 
of 4> in the Hartle-Hawking vacuum by taking the arguments of to the boundary (see 
e.g. |26|,[27| ) .i For example, the boundary Wightman function is obtained by 



G + (x, x') = lim (2vr ){2vr' )Q + (x, r; x', r') (2.10) 

r,r' — >oo ^ ' 

where Q+ and G+ denote the bulk and boundary correlation functions respectively 

G+(x, x') = Tr le- pH 0(x)0(x')] , (2.11) 

g + (x,r;x',r') = (0\<l>{r,x)(l>{r',x')\0) HH . (2.12) 

In the above equations we used the notation x = (t, e) with e denoting a point on S 3 
and the subscript "HH n denotes the Hartle-Hawking vacuum. Going to momentum space 
( |2.10| ) becomes (see Appendix A for a definition of the Fourier transform) 



G+(u,l)= lim (2vr A )(2vr' A )g + (uJ;r,r') , ( 2 .1 3 ) 

r,r — >oo 



where / denotes the angular momentum on S s . The Feynmann (retarded) propagator in 
the bulk leads to the Feynmann (retarded) Green function on the boundary by the same 
procedure. In this paper we will focus on G+ for reasons to be commented on later. 

Since the extended black hole background has two asymptotic boundaries, we can also 
take r and r' to different boundaries. Such points are always space-like separated in the 
bulk and lead to boundary correlation functions of complex time separation (see equation 

G 12 (t) = Tr [e~ m O{t - i/3/2)O(0)] = G+(t - i/3/2) (2.14) 

where we have suppressed the boundary spatial coordinates for notational simplicity. G+ (t) 
is analytic for —(3 < Imt < and the two-sided correlator Gi2(t) can be obtained from 
G + (t) by simple analytic continuationB. 

G+ can be found in terms of solutions to the Laplace equation for <fi following the 
standard free field quantization procedure. Let 

(p = e~ 2Ujt Y I (e)r~ ± 2 1 ip(uJ 1 p; r), 



6 For discussion of boundary Lorentzian correlation functions in the supergravity approxima- 
tion, see also p|MBaM|. 



7 In contrast, the Feynman and retarded functions have singularities in the range — f < Inr£ < 

0. 
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with lj(e) denoting scalar spherical harmonics on S 3 (see appendix A for notations on 
spherical harmonics). The Laplace equation for <p can then be written in terms of the 
tortoise coordinate as a Schrodinger equation for if; 



(-d 2 z +V l (z)-u 2 )ifj = . 



(2.15) 



Vi(z) can be expressed through r(z) as 

V l (z)=f(r) 



j.2 ^ 4r^ 



(2.16) 



For real / > 0, Vi(z) is a monotonically decreasing function of z G (0, oo). Near the 
boundary, 



Vi 



and near the horizon 



Vj oc e 3 







2! — ► + OO 



(2.17) 



(2.18) 



For any given real to the Schrodinger equation (|2.15| ) has a unique normalizable mode 
i/jui, which we will take to be real. We normalize it at the horizon as (6 U is a phase shift) 



z — > +oo . 



As z — > 0, ifj^i has the form 



1> ul KiC{u,l)zt +V + 







(2.19) 



(2.20) 



where the constant C is fixed by the normalization of (|2.19|) . It is easy to check that ip^l 
is even in uj. 

Using the mode expansion of <p m the Hartle-Hawking vacuum, the bulk Wightman 
propagator Q+ ( |2.14|) in momentum space can be written in terms of ip^i as 



g + (uj,l;r, r) = — _ - (rr ) 2 ip ul (r)tp ul (r ) 
which leads to the boundary G+ upon using ( |2.13|) and ( |2.20|) 



(2.21) 



(2v) 



-C 2 (uJ) 



(2.22) 



2u e?" - 1 

G+ is to be evaluated for real u>, I and can be analytically continued to general complex u 
and I. 
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2.3. Analytic properties 

Equations ( |2.21| ) and (|2.22|) indicate that the boundary Yang-Mills theory has a con- 
tinuous spectrum in the large iV and large 't Hooft coupling limit, despite being on a 
compact space. At finite N, the theory should have a discrete spectrum on S 3 . In the 
bulk the continuous spectrum arises due to the presence of the horizon. 

The analytic properties of G+ in the complex t^-plane for a given I can be deduced by 
applying standard techniques of scattering theory to (|2.15| ). The fact that r is a periodic 
function of the tortoise coordinate with a period implies that: [[23] 

1. The poles in the prefactor -^jkrziY) cance l with zeros of C 2 . G+ is analytic at u = 
and uj = neZ. 

2. The only singularities of G+ in the complex w-plane are poles. The locations of the 
poles obey a reflection symmetry: if there is a pole at u>q, then there are poles at 

—Uq, OJq, — Uq. 

These two features are quite generic, applicable to AdS black holes of all dimensions. The 
poles of in the lower half w-plane coincide with those of the retarded propagator Gr 
and correspond to quasi-normal frequencies of the black hole background. Since it is not 
known how to solve the Schrodinger equation (|2.15| ) exactly, the determination of the 
quasi-normal poles is a difficult mathematical problem. 




Fig. 2: Poles for G+(u), I) for I = in the complex w-plane. We use ro = 1, r\ = \/2. 



When / = (or small compared to ^ or w), the problem simplifies and various meth- 



ods [|3^ , |33| , p4| , |35|j3q1 can be used to determine the locations of poles of G + approximately. 
One finds that there are four infinite lines of poles as indicated in fig. [|. The poles in the 
upper right quadrant are given by [53 



27r 47rn 
UJ ~ — V + UJn H — . 

B B 



n 



0,1, 



(2.23) 



The other lines are obtained by reflections, in ( P-23|) is a constant of O(l) (independent 
of v) whose exact value is not relevant here!. 

The quasi-normal frequencies for / 7^ and other dimensions are more complicated to 
find and will be discussed in |24[ . 



3. Boundary Wightman functions and bulk geodesies 

3.1. A semi-classical approximation 



We now develop a "semi-classical" approximation to equation (|2.15|) , in the following 
large v limit 

iv = vu, I + 1 = uk, v 1 , (3-1) 
i.e. we take the mass m of (p to be large and "measure" the frequency ui and angular 



momentum I in units of m. With tp = e uS equation ( |2.15 ) becomes 



-{d z Sf - -d 2 z S + V(z) + \Q(z) = u 2 



(3.2) 



with 



V(z) 



Q(z) = f(r) 



1 1 9/j 

4t-2 4 4^4 



(3.3) 



V(z) 




Fig. 3: The potential V(z) with u = 10 and k = is shown. z c is the turning 
point. Dashed and solid lines indicate the classically forbidden and allowed regions 
respectively. 



The accuracy of equation ( 2.23| ) increases with n. A comparison with the numerical results 



obtained in [32] shows that it appears to work well even for n small. 



8 



With k 2 > the leading order potential V(z) is a monotonically decreasing function 
for z G (0, +00) as indicated in fig. |3|. For scattering sates with u > 0, ( |3.2| ) can be solved 
order by order in 1/^ expansion using the standard WKB method. In the classically 
forbidden region (see fig. |3|), the exponentially decreasing solution can be written as 

^T* 6) W = -i=e^(l + 0(i/- 1 )) (3.4) 

V J K r 

with! 

Z{r) = - [ dr'K r (r'), K r = \^V{r) - u 2 . (3.5) 



/ 

r c in the lower integration limit of (|3.5| ) is the turning point, given by the real positive 
root of the equation 

V(r) = f(r) (l + ^ = u 2 . (3.6) 
For u 2 > 0, equation (|3.6|) has a unique positive root r c > r . Z satisfies the equation 

fZ' 2 ~% + - f u 2 = l (3.7) 

with Z'{r c ) = 0. Note that we have written the above equations in terms of r for con- 
venience. One can equivalently write them in terms of the tortoise coordinate z. The 
expressions in terms of r are more convenient to visualize the analytic continuation to be 
discussed later. 

The expression for ifj^l in the classically allowed region of the potential ( |3.3| ) (i.e. 
for z c < z or tq < r < r c ) follows from the standard connection formula, from which one 
can determine the relative normalization between ifi^^ of ( |3.4j ) and tfj^i of (|2.19| ) to be 

^>fc&) = ^ m > 0, v -> 00 (3.8) 

Ju 



From (|3.8|) we find that in the limit (|3.1j ) the boundary Wightman function can 
be expanded as 

G+(u, I) w 2u e uZ( - u > k) (l + 0(z/ -1 )) + subdominant terms, u > (3.9) 

with 

Z(u,k) = 2 lim (logr- / dr' K r (r')) . (3.10) 



9 The branch cuts for K r on the complex r-plane are chosen so that they do not intersect the 
integration contour in Z. 
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Higher order 1/v corrections in ( |3.9|) can also be obtained from (|3.2| ) using the standard 
WKB procedure. In particular, the term proportional to Q(z) will be important at order 
v . There could also be subdominant terms in (|3.9| ) coming from reflections at other 
(complex) turning points of V{r). 

While equations (|3.9|) -( |3TTC)| ) were obtained for u > and k > 0, they can be analyt- 
ically continued to the full complex u and /c-planes. We will show in section 4 that the 
analytic continuation allows us to probe the region beyond the horizon. 

3.2. Relation with geodesies 

We expect Z(u, k) in ( |3.10|) to have a simple interpretation in terms of bulk geodesies. 
The reason is that in the large mass limit, the propagation of bulk field (p should ap- 
proximately follow geodesic paths. Thus we expect a direct relation between the WKB 
approximation of the last subsection with the geodesic approximation. The scaling in ( |3.1| ) 
simply defines u as the "velocity" in t direction and k as the "angular velocity" on S 3 . 

Due to translational invariance in t and isometries on S 3 , a bulk spacelike geodesic is 
characterized by the integrals of motion 

E = f — , Q = r 2 — (3.11) 

as ds 

where s is the proper distance and 9 denotes the angular coordinate along the geodesic 
motion on S 3 . We treat geodesies which are related by a translation in t and on S 3 as 
equivalent. The geodesic satisfies the equation 

1 fdr\ 2 q 2 1 2 . 



Equation (|3.12[ ) is precisely ( |377| ) with the identification!^ 



(IT 

fZ' = —, u = iE, k = iq. (3.13) 

n r of ( |3.5| ) can be identified as the proper velocity of the geodesic along the r direction. 
Thus Z(u, k) can be associated with a (complex) spacelike geodesic with constants of 



10 The sign choice for the first expression below corresponds to having the geodesic moving 
away from the boundary, while those for the last two equations are for convenience. The fact i 
in relating the boundary "velocities" u, k to the velocities E, q of the bulk geodesic is due to that 
the geodesic is spacelike. 
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motion E = —iu and q = —ik, which starts and ends at r = +oo.0 More explicitly, one 
finds that Z(u, k) can be written as 

Z(u, k) = —Et(E, q) - L(E, q) + qd(E, q) (3.14) 

where L(E, q) is the (regularized) proper distance of the geodesic, t(E, q) is the time 
separation and d(E, q) is the proper distance on S 3 between the final and initial points, 

L{E, q) = 2 lim f / £ = - log r 

f°° dr 

t{E, q) = 2E = (3.15) 



d(E,q) = 2q 
Also note the relation 



dr 



r 2 f + E 2_f q 2 



87 r)7 
m = -Wq) -^ = d{E,q) (3.16) 

which shows that L(t, d) and Z(E, q) are related by a Legendre transform. 

Note that E and q do not specify a geodesic uniquely. ( |3.15| ) defines a complex geodesic 
with a choice of root r c (E,q) of equation ( |3.6|) as the turning point and a contour from 
r c (E,q) to +oo. For the same value of E, q, a different choice of the root or a different 
contour which cannot be smoothly deformed into the previous one defines a different com- 
plex geodesic. The identification ( |3.14j) with the boundary Z(u, k) selects a specific one 
among them. 

In the above discussions we have concentrated on the Wightman functions^!. A sim- 
ilar relation with bulk geodesies can be established for retarded and Feynmann functions 
in momentum space, whose story is more complicated since their dependence on v is not 
uniform. Special care is needed when v is an integer^!. This makes the large v limit more 
subtle. Even at the leading order, one has to take into account of an infinite number of 
classical paths in the WKB approximation [ 37)] . Fortunately, all these additional compli- 



cations arise due to the asymptotic behavior of / ~ r 2 near the boundary of spacetime 
and do not seem to give additional insight into the question of physics beyond the horizon. 



11 Note that depending on the values of E and q the starting and end points can be on the 
same or different boundaries. 

12 Wightman functions already contain all the information of the theory. For example, Feyn- 
mann and retarded functions can be obtained from them. 

13 This can be seen also from zero temperature correlation functions. 
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3. 3. Coordinate space correlation functions 

We now look at the Fourier transform of G+(u,l) to the coordinate space correlator 
(see Appendix A for notations) 

i 00 r°° a 

G + (t;e,e') = i ^^2(/ + l)G(e-e / ) J —e**G+{u,l). (3.17) 



The two-sided correlator (|2.14|) can be obtained from (|3.17| ) by taking t — > t — i^, while 
the Euclidean correlator G^(r; e, e') can be obtained by taking t = —ir with < r < (3. 

In the large v limit (|3.1|) , using (|3.9|) we can approximate the sum over I in ( [3.171) by 
an integral over k 

,,3 f-oo 

G+(t,9)fa —= -/ dudkke- lvut+lvke 2ve vZ{u > k) (3.18) 

9>-K 6 is\n9 J_ 00 

where 9 = cos _1 (e-e / ) and we have extended the integration range for k to (— oo, oo) using 
that Z is an even function of k. ( |3.18j ) can be evaluated by the method of steepest descent 
with the saddle points determined by 

Using equations ( |3.16|) and ( |3.13|) we find that ( |3. 19|) become 

t = t(E,q), 9 = d(E,q) (3.20) 

i.e. bulk geodesies with end point separation given by (t, 9) appear as saddle points of 
( |3.18|) . Since from (|3.14|) the regularized geodesic distance L(t, 9) and Z(u, k) are related 
by a Legendre transformation, one finds that 

G + (t, 9)^J2 2 » J * (^) 2 e ~ VU i 1 + (3-21) 

i 

where i sums over the saddles along the steepest descent contour. The Jacobian J is due 
to the Gaussian integration around the saddle points. 
Some remarks: 

1. J can be interpreted as the density of the geodesies. It is proportional to the Van 
Vleck-Morette determinant for the geodesies. One can check that ( |3.21| ) agrees pre- 
cisely (including the prefactor J) with the expression obtained directly from the 
geodesic approximation to the coordinate space path integral 

G(z,r;z',r')= £ e**" 3 , (3 22) 

paths 
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after taking the end points to the boundary. 

2. In the standard geodesic approximation to ( |3.22| ), it is often a subtle question in 
Lorentzian signature to determine which geodesies contribute to the sum ( |3.21|) . In 
our approach, the steepest descent approximation of the Fourier transform ( |3.17| ) gives 
a precise prescription for determining the sum. 

3. Higher order terms in (|3.21| ) can be computed systematically in our approach. 



4. Black hole singularites in Yang-Mills theory 

In this section we discuss how information about the black hole singularity can be 
extracted from Z(u, k) using its relations with the bulk geodesies developed in the last 
section. For simplicity, we restrict our discussion to k = 0, in which case the corresponding 
bulk geodesic has zero angular momentum. The discussion for general k is more involved 
and will appear in |[23|| . The k = case already captures many of the essential elements. 



4-1. Probing the physics beyond the horizon and near the singularity 

We first consider the analytic continuation of Z(u, k = 0)0 to general complex u. We 
will see that the analytic continuation allows us to probe the regions beyond the horizon 
and near the singularity. 



As discussed in section 2.3, when u > 0, the turning point r c (u) in (|3.10| ) lies outside 
the horizon, i.e. r c > tq. The integration contour runs along the positive real r-axis from 
r c to +oo. The analytic continuation of Z(u) to the full complex it-plane involves the 
following two aspects: 

1. Analytically continue the turning point r c (u) from that for real u > 0; 

2. Smooth deformation of the integration contour as r c moves on the complex r-plane. 
Both steps have some subtleties, which we now discuss in detail. 



14 In the following, we will simply write it as Z(u) and similarly use G+(cu) for G+(ui, 1 = 0). 
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S_ 

Fig. 4: The structure of branch cuts of Z(u) and r c (u) for tq = 1, n = 2. At 
finite v, the branch cuts become the pole lines of G+(uj,1 = 0) as in fig. |||. The 
asymptotic regions are labelled by S± or B±, indicating whether the turning point 
for the corresponding u approaches the singularity (<S±) or the boundary (B±). 



The analytic continuation of r c from the u > region is not unique, since r c (u) has 
branch points in the complex it-plane at which it merges with other solutions of ( |3.6| ). 
These are also branch points of Z(u). When k = ( |3.6| ) is a quadratic equation for r 2 and 
r c (u) coincides with other roots when 

(u 2 - l) 2 + 4^i = . (4.1) 



From (|4.1| ) we find that r c (u) has four branch points at 

2tt 2tt s 
u = ±(r 1 ±ir )=±=-,±— . (4.2) 

For r c and Z to be single-valued on the it-plane, branch cuts have to be specified. Dif- 
ferent choices of the branch cuts correspond to different ways of performing the analytic 
continuation. The locations of the branch cuts cannot be determined from the integrals 
( p.l0| ) or ( |3.15| ) alone. To determine them we need to use analytic properties of G+(u). 
As discussed around equation Q2.23D , the only singularities of G+(u) at finite v are four 
lines of poles located at 

27T Wn 47T71 , , . 

+^+ n = 0,l,..-. 4.3 

B v vB 



and the reflections of (|4.3|) with respect to the real and imaginary u axes. In the large 
v limit, since the spacings between poles go to zero, these lines of poles become branch 
cuts of Z(u)B This determines the directions of the branch cuts to be along the radial 
direction from each branch point to infinity (see fig. HI). 



15 The simplest example exhibiting this behavior is a Gamma- function Y{yz). The function has 
poles along the negative real axis for finite v. In the large v limit, upon using the Stirling formula, 
the pole line is replaced by a branch cut. 
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Fig. 5: A radial spacelike geodesic can be described by a particle of energy —u 2 
moving in the potential U = —f. The horizon is at r = ro- For u 2 < 0, the turning 
point lies inside the horizon. 



With the branch cuts precisely specified, r c {u) can now be uniquely determined from 
that for u > 0. In particular, fig. [| implies that the analytic continuation should be done 
through the region around u = 0. To to be definite, let us concentrate on real u 2 . In 
this case, a convenient way to visualize how the turning point changes with u is to treat 
equation ( |3.12|) as the motion of a one-dimensional particle of energy E 2 = —u 2 , moving 
in a potential^ 

U = —V = -f , (4.4) 

as in fig. |5|. For real u, r c > ro, i.e. the turning point lies outside the horizon, while for 
u pure imaginary, r c < ro and the the turning point lies inside the horizon. One can also 
solve (|3.6|) explicitly and r c (u) is given by the positive branch of the equation 



r? = ,/„ + (i^i) 2 -i^I. (4.5) 



With the turning point specified, one can now find the bulk geodesies corresponding 
to various values of u from simple geometric considerations. For example, real values of u 
correspond to radial geodesies in the Euclidean section of the spacetime (see fig. ^), since 
from (|3.13|) and ( ft. 11 ) t is pure imaginary along the geodesies. Pure imaginary values of 



u = iE (E real) correspond to real geodesies in the Lorentzian section of the spacetime, 
which connect two asymptotic boundaries (see fig. 0). The turning point of such a geodesic 
lies in region II for E > and in region IV for E < (see also fig. |l|). 



16 Note the potential V is inverted since we work in the classically forbidden region of the 
Schrodinger problem (|2.15| ). 
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Fig. 6: Radial geodesies in the Euclidean section of the spacetime which corre- 
sponds to real values of u. The Euclidean section of the r — t plane is a disk with 
it as the angular coordinate and the origin of the disk at r = ro- The solid circle 
is the boundary. Geodesies with u > (i and ii) and u < (iii) are schematically 
plotted. Geodesic ii correspond to the large u limit, in which case the tuning point 
is close to the boundary and the end points of the geodesic are nearly coincidental. 




Fig. 7: Radial geodesies in the Lorentzian section of the spacetime corresponding 
to pure imaginary values of u = iE are schematically plotted. Geodesies i and ii 
have E > while iii has E < 0. Geodesic ii correspond to the limit E — > +oo, in 
which case the tuning point is close to the singularity and the geodesic becomes 
nearly null. 

The dependence of r c on u illustrates some interesting features in the relation between 
bulk and boundary scales. For real u — > ±oo, the turning point is given by 

r c w \u\ — > +oo (4-6) 

i.e. the turning point approaches the boundary. In this limit, the end points of the geodesic 
becomes nearly coincidental. When u decreases, r c also decreases. The turning point r c 
reaches the horizon for u = 0. This behavior reflects a familiar feature of the AdS/CFT 
correspondence, called IR/UV connection [p]JS^|, which relates long distances in the AdS 
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spacetime to high energies in the boundary theory. The turning point r c moves inside the 
horizon when u moves along the imaginary axis from the origin. Let u = iE. Then as \E\ 



i.e. the turning point approaches the singularity. Thus when dealing with physics inside 
the horizon, there appears to be a new feature. To probe deeper inside the horizon requires 
larger E. Since the singularity may heuristically be considered as the UV of the bulk, we 
find a UV/UV connection. It is important to keep in mind that inside the horizon, r plays 
the role of the time coordinate. This indicates that the "time" inside the horizon is indeed 
holographically generated from the boundary Yang-Mills theory. 

The above discussions can be easily generalized to all complex values of u using equa- 
tion ( ^.6| ). In particular, from ([OJ), \r c \ ^ requires that \u\ — > oo, due to the fact that 
/ blows up at the singularity (large curvature effect). Conversely, \u\ — > oo implies either 
\r c \ ^ or \r c \ — > +00. Thus along different directions to infinity in the complex w-plane, 
the turning point either approaches the boundary or the singularity. The branch cuts 
in fig. |] divide the complex infinity of the it-plane into various asymptotic regions. The 
regions which correspond to the singularity or the boundary are indicated in fig. £|. Near 
the real u axis, the turning point approaches the boundary as \u\ — > +00. As \u\ — > +00 
near the imaginary u axis, the turning point approaches the singularity. 



increases, r c decreases (see equation (P~5|)). For \E\ — > +00, we find that 



\E\ 







(4.7) 





Fig. 8: The integration contours in the complex r-plane for (a): u > 0, (b): u < 0, 
(c): u = iE, E > 0, (d): u = iE, E < 0. The solid dot indicates the pole r = ro of 
the integrand (horizon). 
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We now look at the second aspect of the analytic continuation of Z(u). This involves 
the deformation of the integration contour in (|3.10| ) (or (|3.15| )) as r c (u) moves in the com- 
plex r-plane. As the contour is deformed, one needs to be careful about the contribution 
of the pole at r = tq of the 1/ / factor in the integrand. To find the contribution of the 
pole, it is enough to consider when the turning point is close to the pole0. This happens 
when \u\ is small, in which case a prescription for the contour deformation can be obtained 
from the requirement that Z(u) be analytic at u = 0. The contours for other values of u 
can then be obtained by continuous deformation without further subtlety. The resulting 
contours for real and pure imaginary u's are plotted in fig. |8|, from which the contribution 
of the pole can be readily obtained. 

For real u < 0, the contribution from the pole to t(u) of ( |3.15| ) is —i/3. This is precisely 
the period of the complex time and matches with the geometric picture of fig. [| When 
u = iE for real E, the imaginary part of t(u) solely comes from the pole contribution and 
we find 

Tmt(iE) = - — , E real . (4.8) 

This shows that the end point of the corresponding geodesic lies in the other asymptotic 
boundary, i.e. region III of fig. [I], again consistent with fig. [F[ 

To summarize, through the function r c (it), we establish a correspondence between the 
complex u- and r-planes. G+(vu) evaluated at u in the large v limit probes the black 
hole geometry near r c (u). The boundary correlation functions encode not only the bulk 
geometry outside the horizon, but also regions beyond the horizon and near the singularity. 

We emphasize that the proper identification of the branch cuts for Z(u) and r c (u) 
(which in turn depends on the knowledge of poles for G+(u)) is crucial for our conclusion 
above. Different choices of the branch cuts could lead to completely different physical 
pictures. For example, a different analytic continuation procedure may lead to a r c (iE) 
that for real E is not given by the positive branch of (|4.5| ), but some other root of the 
turning point equation ( |3.6[ ). In that case one cannot associate the geodesies in fig. [7| with 
G+(u) and it is not clear one could probe the physics beyond the horizon. 

Given Z(u) from the large v limit of the boundary G+(yu), the problem of finding the 
bulk metric is essentially a classical inverse scattering problem. Due to the large number 
of isometries of the background, the problem is effectively one-dimensional, being that of 

17 / also has other poles at — ro and Uri. The analytic continuation procedure is such that 
the turning point never coincides with these poles. 
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a particle moving in (O). The equivalent one-dimensional problem can be phrased as 
follows. Consider sending a particle toward the potential from r = oo and waiting for it to 
come back. L(E) is then the (regularized) time interval for this scattering process. With 
the knowledge of L(E) for all values of E, one can in principle reconstruct the potential 
( |4.4| ). At a given E, L(E) probes the behavior of the potential ( f4.4| ) near the turning point 
r c (u). 



4-2. Manifestations of singularities in boundary theories 

We have found that the geometry around the black hole singularity is encoded in the 
behavior of G+(iv) near the imaginary infinity. We now examine the manifestations of the 
singularity explicitly. 

The integrals in (|3.15|) can be evaluated explicitly and one 
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In ( |4.9| ) the branch cuts of the logarithms are chosen to be straight lines extending radially 



from ±^ and ±4£ to 00, as follows from the discussion in previous subsection (see fig. |j). 



Expanding ( ^.9| ) for large \u\ we find that L(u) and t(-u) of (|3.15|) can be written as 
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The expressions below were obtained before in [12| in a somewhat different form. 
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and to an d [u] are given by 
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(4.13) 



where B± and S± denote asymptotic regions in fig. ^ whose corresponding turning point 
approaches the boundary and the singularity respectively. The values of to fo r various limits 
can be easily understood from the geometric pictures of the geodesies in fig. || and fig. [7[ 
When u — > ±oo, the end points of the Euclidean geodesies become nearly coincidental. 
The value —i(5 for u — > — oo is precisely the full period of the Euclidean circle. When 
u — > ±zoo, the bouncing geodesies in fig. ^ become nearly null and the corresponding 
values of to i n Q4. 13|) are twice of those in ( |2.5|) . 

Equations (|4.12| ) implies that as uj = vu — >• ±zoo, the boundary correlation function 
behaves as 



G + (uJ = 0) 



1 



t(TM) ! 



.UJ 



1v i w [±4_M 



l + O 



1 



(4.14) 



where the upper (lower) sign corresponds to uj — > +zoo (—zoo). Note that the correlation 
function decays exponentially along these directions. For uj — > ±oo near the real axis, we 
find 



(ir(i+o(^)) 



G+(u,l = 0) 



F(T>)F 
l 



UJ 



(-1) 



2// 



+oo 



-oo 



(4.15) 



W ru e/3w C 1 + (^)) 

While equations ( 4.14 ) and ( [4.15|) were derived in the large v limit, they should hold 
for finite u, since the \u\ —>■ oo limit should coincide with the limit \uj\ = v\u\ — > oo 
regardless of the value of v. Note that ( fll5| ) is precisely what one would expect of the large 
frequency behavior of a conformal field theory at finite temperature^!. The exponential 
falloff in ( |4.14j ) reflects the presence of a curvature singularity in the bulk. The falloff is 
controlled by the complex parameter B (introduced in (|2.6| )) which characterizes the black 
hole geometry. 
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LO 



+oo, one expects the correlation function to recover the zero temperature 



result. The second line of equation (4.15) follows from the general properties of the Wightman 
function at finite temperature. 
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4-3. Generalizations to nonzero angular momentum 

The above discussions can be extended to Wightman functions of nonzero angular 
momenta. We summarize some main results here, leaving detailed discussions to [ |23| : 

1. For boundary angular velocity k real, which corresponds to geodesies of pure imaginary 
angular momentum, the structure of the branch cuts for Z(u, k) is similar to fig. f|. 
The locations of the branch points and the directions of the branch cuts depend 
nontrivially on k. At finite u, the branch cuts become lines of poles of G+(u, I). 

2. For bulk geodesies with real angular momentum q, an important new feature appears: 
there exist geodesic orbits with constant real r. The existence of such orbits leads 
to the appearance of virtual states (if there is an orbit lying inside the horizon) or 
bound states (if there is an orbit lying outside the horizon) in the Schrodinger problem 
Q2.15Q . These virtual states or bound states lead to two new lines of poles of G+(u>, I) 
along the imaginary w-axis for pure imaginary I. Thus Z(u, k) has two new branch 
cuts along the imaginary u-axis for k = iq pure imaginary. 

3. It remains true that the turning point approaches the boundary (singularity) for \u\ —> 
oo along the real (imaginary) axis. Furthermore, for any fixed /, in the large u limit, 
equations ( |4.14j )-( fi.l5| ) remain valid. 

4. In the limit that q goes to zero, the branch cuts for Z(u, iq) along the imaginary axis 
move to infinity and fig. ^| is recovered. The fact that there are branch points at 
u = ±zoo at q = leads to interesting behaviors in the expansion of Z(u, k) around 
k = 0. For example, let u = iE with real E, then one finds that to leading order in 
the limit E — ► +oo, Z(u, k) has the following small k expansion^ 

Z(iE, k) « + 2 log | + 2 + -L f; ai (k 2 E 2 ) 1 + ■■■ (4.16) 

i=i 

with 

7T 37T 1 

Oi = -— — , a 2 = — — , a t 5^ (4.17) 

2/i2 16//2 

where \x dependence can be deduced based on dimensional analysis. Note that the 
expansion parameter for small k is given by k 2 E 2 and the derivatives over k at k = 
become divergent in the large E limit. We will see in the next section that Q4.16| ) 
leads to divergences in certain gauge invariant observables in the boundary theory. 



20 Note that the equation below is obtained by first doing a small k expansion and then taking 
leading order terms in 1/E. 
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We also mention by passing a few other generalizations: 
5. The discussions of this section can also be generalized to an AdS^+i black hole of 
dimension d > 2 [f23J . All the essential features for AdSs black holes carry over to other 
dimension d > 3. A special case is a BTZ black hole in AdS3, which has an orbifold 
singularity and the story is somewhat different. In the BTZ case, the corresponding 
Schrodinger equation fl2.15|) can be exactly solved and the relation between the large 



v limit of Wightman functions and bulk geodesies can be explicitly verified. 
6. The subdominant contributions in ( |3.9| ) can also be worked out using a more sophis- 
ticated WKB method involving more than one turning point. One can show that an 
infinite number of subdominant contributions become important at the branch cuts0 



where they add up to produce the poles of G+(u,l). In |24j we use this property to 
derive the positions of poles of G+ for general / and dimension in the large v limit. 
7. While we have not examined it in detail, it is interesting to compute the higher order 
1/u corrections in ( |3.9| ). In particular, the function Q(z) ( |3.2| ) will start contributing 
at the order 0(l/u). Since Q(z) becomes singular at r — > 0, it would be interesting 
to see whether it yields new manifestations of the singularity in the boundary theory 
correlation functions. 

4-4- Coordinate space correlators and alternative indications of curvature singularities 
In this subsection we consider the Fourier transform of G+(c<;, I) to coordinate space. 



To make connection to the result of 12 , we consider 



G 12 (t, 6 = Q) = G + [t- % ^6 = ^j (4.18) 

which can be obtained from ( |3.17| ) by taking t — > t — ^ and corresponds to inserting 
operators on two different boundaries. We restrict to 9 = for simplicity. In the large 
v limit, G\2 can be evaluated in exact parallel of the discussion of section 3.3. One can 
approximate the sum over I by an integral 

/ duAhh 2 p-ivut-^vuP uZ{u,k) 



G 12 {t)tt—\ dudkk 1 e -^t-^up 2l/e vz (u ,K) ( 419 ) 



and perform the integrals using the saddle point method. 



21 In other words, the branch cuts are anti-Stokes lines for an infinite number of subdominant 
exponentials. The Stokes lines can be obtained from the branch cuts by rotating ±-|. A simple 
example which also exhibits this phenomenon is the Gamma- function (see e.g. [f39f ). 
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Fig. 9: The contour plot for the real part of Z(u,k = 0) — iu(t — i^) for i < § 
in the complex u-plane. There is a saddle on the imaginary axis and two complex 
ones. The steepest descent contour is also shown in figure. The contour does not 
pass through the saddle on the imaginary axis even if it dominates. 

From the discussion of section 3.3, the saddles of (|4.19| ) correspond to geodesies whose 
end points lie on different boundaries and are separated in time by an amount t and no 
separation on S 3 . The saddle for the /c-integral is simply k = 0. The saddles for the 
it-integral are given by 

8Z(uk = 0) = p 

du 2 y ' 

Note that Z(u, k = 0) can be obtained from (f4~9p. The solutions to (fO0|) can be visualized 
conveniently on the contour plot of the real part of Z(u, k = 0) — iu(t — i§) * n ^ ne complex 
•u-plane. See fig. |^. In the figure we also indicate the steepest descent paths to which the 
integration contour of ( f4.19|) can be deformed. 

The dependence of the saddle point structure on t can be summarized as follows. For 
t < t c = |, there are three saddles, as indicated in fig. |^. The one on the imaginary axis 
corresponds to a real geodesic in Lorentzian black hole spacetime with a turning point 
inside the horizon (see fig. 0). This is the bouncing geodesic discussed by [^2| . We will 
refer to this saddle as the bouncing saddle below. The other two saddles describe complex 
geodesies, which do not seem to probe the physics beyond the horizon. As t approaches 
t c , the bouncing saddle moves to infinity along the imaginary axis and the turning point 
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of the geodesic approaches the singularity. For t > t c , the bouncing saddle disappears^. 

From the steepest descent contour, we conclude that the bouncing geodesic does not 
contribute to coordinate space correlation functions. This result was obtained in [|12] by 
analytical continuation from Euclidean signature!!. Here we confirm their result. 



4-5. New observables of the boundary theory and manifestations of the curvature singularity 



While the bouncing geodesic does not contribute to ( [4. 18| ) directly, from momentum 
space correlation functions, we can easily construct observables in the boundary theory 
which are approximated by the bouncing geodesies in the large v limit. 



C 2 











Fig. 10: The integration contours for Gi2(t) (Ci) and H\2{t) (C2) 



Let us start with two-sided correlator Q4.18| ) with coincidental spatial coordinates 

G 12 {t) = Tr [e-P H (D(t - i/3/2, e)O(0, e)] 

du _ iu „„ , , N (4.21) 



where the contour C\ is along the real w-axis and 



1=0 JLl 



G 12 (u,l) = e-^G+(u,l) . (4.22) 



22 For t > t c , the two complex saddles wind around the branch points at ^ and — = respectively 
and approach them as t — > 00. 

23 In [p.2| it was argued that the information regarding the bouncing geodesies and thus the 
singularity can nevertheless be obtained by analytic continuation in the large v limit, since there 
exist certain values of t for which the bouncing saddle merges with other saddles. 
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From ( 4.14|) -( f4.15| ) and the comments in sec. 4.3 regarding their generalizations to nonzero 
/, we find that the large u behaviors of G^w, I) are given by 



(4.23) 



We now construct new observables in the boundary theory 

duo _„•„,. 



1 /" 

z=o JC 



2tt 



G12W) 



(4.24) 



where the contour C2 are along the imaginary axis. H\2(t) can be defined for — ^ < r < ^ 



due to the exponential falloff (]4.23|) of Gi 2 (^) along the imaginary axis. Note that H\2(r) 
is gauge invariant by definition. While its existence depends on the asymptotic behavior of 
Gx 2 (u), I) along the imaginary w-axis, it is an object which can in principle be intrinsically 
defined in the boundary theory!^. 

In the large v limit (|4.24|) can be evaluated in exact parallel as (^.18| )-( ^.19|) by 
approximating the sum over / by an integral and performing the saddle point approximation 
in the resulting integrals. The only difference is that the integration contour for uj is now 
along the imaginary axis rather than the real axis. Thus instead of picking up the two 
complex saddles in fig. [| (|4.24|) is given by expanding around the bouncing saddle on the 
imaginary axis. 

As 5t = t c — t — ► 0, the bouncing saddle moves to the imaginary infinity and the 
turning point r c of the corresponding geodesic approaches the singularity. In this limit, 
we can use ( f4.16|) to approximate Z(u,k). Then the integrals for Hi 2 can be written as 

(u = iE) 



H 12 (r) 




,ESt + 2vlog- + -^J2 a i ( eE2 ) 1 



1=1 



(4.25) 



24 One can of course define Hi2(r;e,e') for non-coincidental points on S 3 . We do not expect 
divergences for them. They will be discussed in [ p3| . 
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where the saddle for the /c-integral is k = and for the E integral is 

E c « j t -> oo . (4.26) 

Note that in ( f4.25|) we have only worked out the qualitative 5t dependence in the small 
St limit rather than attempting a precise evaluation. The order of limits should be first 
v — > oo and then St — > 0. In evaluating the high order 1/u terms in ( [4.25| ), there are two 
other sources of 1/v corrections that we did not take into account: 1/u corrections in ( |3.9| ) 
and those in turning the sum over / in ( 4.24| ) to an integral. Since we are only interested 



in power counting, apart from magic cancellations, this will not affect the qualitative 
behaviors in ( f4.25|) . The divergences in l/v n terms arise from the /c-integral and are due 
to the structure of the small k expansion in ( f4.16|) . 

Note that (|4.25|) is precisely what one expects of a bouncing geodesic as its turning 



point approaches the singularity. From the bulk point of view, one expects the contribution 
from the geodesic should take the formil 



x ' R 



1 + ( 4 - 27 > 



v 

n=l 



where L is the regularized the geodesic distance, J is the Van Vleck-Morette determinant. 
Higher order terms in ( fj.27p arise from cubic and higher order terms in the sum over paths 
around the geodesic and R n can be expressed in terms of bulk geometric quantities in 
the form of components of curvature tensors (and their derivatives) integrated along the 
geodesic. In general the explicit expressions for R n are very complicated (see e.g. [|IJ). On 
dimensional ground, one expects that in the limit that the turning point approaches the 
singularity R n ~ — , with e the proper time between the turning point and the singularity. 
This leads to 



3C 



jr^ji^ll + E I . (4.28) 



71=1 



Equation ( 4.28[) precisely agrees with ( |4.25| ) since from the metric fl2.2f) 



where we have used ( |4.7|) and (OB). 



25 Similar argument was also used in JT 
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To summarize we have constructed new gauge invariant observables ( 4.24 ) in the 
boundary theory that are sensitive to the physics beyond the horizon and in particular 
their behaviors near r — > | precisely reflect the curvature divergence of the singularity 
These observables are related nonlocally in time with ( [4.1<j| ). It would be interesting to 
better understand their meaning in Yang-Mills theoriesi!. 



5. Discussions: Resolution of black hole singularities at finite N ? 

In this paper we established a direct relation between space-like geodesies in the bulk 
and the large operator dimension limit of the boundary Wightman functions G+(u, I) in 
momentum space. The results present an intriguing picture on how physics beyond the 
horizon is encoded in the boundary theory. In particular, it gives a clear indication that 
the "time" inside the horizon is holographically generated from the thermal Yang-Mills 
theory. 

The poles of G+(u, I) separate the asymptotic region of the complex w-plane into 
several sectors (see fig. ^] and fig. |]). The sectors near the real axis describe the physics 
near the boundary while the sectors near the imaginary axis describe the physics near the 
singularity^. We found the following signals of the singularity in the boundary theory: 

1. G+(u;,Z) falls off exponentially for u — > ±ioo (see equation ( [1141 ) or ( |123|) ). The 
falloff is controlled by the complex parameter B ( |2.6| ), which characterizes the black 
hole geometry. 

2. We constructed new observables #12 (t) (equation ( [4.24] )) in the boundary theory 
which are related nonlocally in time with coordinate space Wightman functions. The 
curvature divergence of the singularity is reflected in the divergences of Hi 2 (t) as 
t — > While the leading order divergence of #12(7") (i.e. the prefactor in ( |4.25|) ) 
can be attributed to 1. above, the divergences in higher order 1/u terms are due to 
more delicate behavior of Z(u, k) for small k and u — > ±ioo. 



26 It seems possible to define a new theory whose momentum space correlation functions evalu- 
ated at it) are given by G+(iuj). #12(1") would correspond to "Euclidean" Green functions in this 
new theory. 

We note that this conclusion does not depend sensitively on whether the curvature singularity 
lies in the Lorentzian section of spacetime. For example, if one regularizes / in (|2.2|) by / = 



r + 1 — r^X^ w ith a small e to move to the singularity to the complex r-plane [15], the picture 



would remain exactly the same, provided this does not change the pattern of quasinormal poles 
significantly, which should be the case for e small. 
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While the above results were derived in the large v expansion, the essential features should 
persist for finite v. For example, we expect equations (|4.14j) and (f4.16| ) should be valid for 
finite v as wellS 

We now comment on the implications of our results in resolving the black hole singu- 
larity. 

The rich analytic behavior observed for G+ in the complex w-plane is tied to the fact 
that in the large N limit, the boundary theory has a continuous spectrum, even though 
it lives on a compact space. In the bulk, the continuous spectrum arises because of the 
presence of the horizon. In Yang-Mills theory, the continuous spectrum should be related 
to the fact that in the high temperature phase, typical states in the thermal ensemble have 
energy of order Y 2 .H 

At finite N, no matter how large, the boundary theory on S 3 has a discrete spectrum. 
In particular, the finite temperature Wightman function should have the form 

G+(u) =27rY,e- 0Em p m n5(u-E n + E m ) (5.1) 

m,n 

which is a sum of delta functions along the real w-axis, where m, n sum over the physical 
states of the theory. G+(u) in equation ( |OD does not have an unambiguous continuation 
off the real axis. In particular, the procedures of analytically continuing G+ to complex 
uj and taking the large iV limit do not commute. Equation ( f4.14|) arises by taking the 
large iV limit first and then doing the analytic continuation. This appears to imply that at 
finite N, geometric notions associated with a black hole, such as the event horizon and the 
singularity, no longer exist. This is not surprising since the black hole geometry arises as a 
saddle point in the path integral of the boundary theory in a 1/N expansion. If one does 
not use such an expansion, the geometric notions lose their meaning. Thus the singularity 
appears to be resolved at finite N. 

The above arguments, however, do not tell us how the singularity is resolved. There 
are several possibilities according to which the singularity can be resolved in AdS/CFT: 

28 This can be checked explicitly using other approximations. 

29 These states have degeneracies of order e cN for some constant c in free theory. When turning 
on interactions, one expects the degeneracy is lifted and the energy levels have spacings of order 
£ -cN This gives rise to a continuous spectrum in the infinite ./V limit. We thank O. Aharony, 
M. Douglas and S. Minwalla for a discussion of this point. Also note that in the low temperature 
phase, the Yang-Mills theory has a discrete spectrum even at N = oo. 
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I. The singularity is already resolved by a'-effects in perturbative string theory, i.e. at 

finite 't Hooft coupling and infinite N. 
Ha. The singularity is resolved only at finite N. But at infinite iV there is a large iV phase 

transition at certain value of the 't Hooft coupling, 
lib. The singularity is resolved only at finite iV and there is no large iV phase transition 

for any 't Hooft coupling. 
To see which possibility is realized, it is important to investigate whether the signals of 



the singularity found this paper and in |T2[ persist to weak coupling in the large iV limit. 
If the answer is yes, it would strongly suggest possibility lib above. This would be a very 
desirable situation since one would then be able to study the black hole singularity in 
string theory using weakly coupled Yang-Mills theory and focusing on the large N limit. 
If the answer is no, then both I and Ha are possible. In the event that Ha is realized, one 
should still be able to detect signals of the singularity at weak coupling, even though the 
precise signals may not be directly obtainable from the results at strong coupling. 

In any case, we believe the results in the paper should provide a valuable guide for 
understanding the black hole singularity in AdS/CFT. 

Finally, it would be interesting to apply the techniques we developed in this paper to 
other backgrounds, like charged or rotating black holes. 
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Appendix A. Fourier transform on S 3 

A complete set of scalar harmonics on S 3 can be written as Yi mm >(e) transforming 
under (1/2, 1/2) representations of 50(4) = SU(2) x SU(2) with -1/2 <m,m' < 1/2. We 

29 



use e to denote a point on S 3 and I = (/, m, m') to denote the full set of indices. Yj is 
normalized so that 



with 



Also note that 



/ Y I {e)Y J {e) = 5 I j 

E ^(ei)^(e 2 ) = ^2(1 + l)C,( ei • e 2 ) 



sin(/ + 1)0 

CHcos0) = — 

smy 



v| 3 y J = -/(/ + 2)y 7 



where V| 3 is the Laplace operator on ,S 3 . 

Consider a correlation function G(e, e') on ,S 3 which only depends on the geodesic 
distance between two points e and e' . It can be expanded as 



1 00 



1=0 



= Y,G{l)Y I {e)Y\e>). 
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